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ABSTRACT 

Supermassive black hole binaries (SMBHBs) are products of galaxy merg- 
ers, and are important in testing A cold dark matter cosmology and locat- 
ing gravitational-wave-radiation sources. A unique electromagnetic signature of 
SMBHBs in galactic nuclei is essential in identifying the binaries in observations 
from the IR band through optical to X-ray. Recently, the flares in optical, UV, 
and X-ray caused by supermassive black holes (SMBHs) tidally disrupting nearby 
stars have been successfully used to observationally probe single SMBHs in nor- 
mal galaxies. In this Letter, we investigate the accretion of the gaseous debris 
of a tidally disrupted star by a SMBHB. Using both stability analysis of three- 
body systems and numerical scattering experiments, we show that the accretion 
of stellar debris gas, which initially decays with time oc t~ 5 ^ 3 , would stop at a 
time T tr ~ rjT^,. Here, t] ~ 0.25 and 21 is the orbital period of the SMBHB. 
After a period of interruption, the accretion recurs discretely at time T r ~ £TJ,, 
where £ ~ 1. Both rj and £ sensitively depend on the orbital parameters of the 
tidally disrupted star at the tidal radius and the orbit eccentricity of SMBHB. 
The interrupted accretion of the stellar debris gas gives rise to an interrupted 
tidal flare, which could be used to identify SMBHBs in non-active galaxies in the 
upcoming transient surveys. 
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Introduction 



Supermassive black hole binaries (SMBHBs) are pre dicted by the hierarchical galaxy for- 
matio n model in A cold dark matter (ACDM) cosmology ( Begelman et alll980 ; Volonteri et al 
20031 ). After a SMBHB at the center of merging systems become hard, it would stall at the 
hard radius for a timescale even lo nger than the Hubble time if spherical two-body relaxation 
dominates (IBegelman et al.lll980l ). However, recent investigations suggested that the hard- 
ening rates of SMBHBs can be boosted and SMBHBs may coalesce within a Hubble time 
either due to various stellar dynamical processes other than the spherical two-body relaxation 



(Yu 



2002 : IChatteriee et al.ll2003l: Merritt fc Poon 



20081 ). or due to gas dynamics (IGould fc Rix 



2004 



2000 



Berczik et al. 



Liu et al. 



200J 



2006: 



& IColpi 



Sesana et al. 



i & Dotti 



2009 



and references therein). 

The strong gravitational wave (GW) radiation generated by coalescing supermassive 
black holes (SMBHs) is the main target of the GW detector the Laser Interferometer Space 
Antenna (LISA) and of the Pulsar Timing Array (PTA) program. Because of the poor 
accuracy of both LISA and PTA in locating GW radiation sources, it is of key importance to 
detect electromagnetic counterparts (EMCs) of GW radiation sources. Identifying SMBHBs 
by their EMCs is also essential to constraining the poorly understood galaxy-merger history. 
Several EMCs have been suggested in the literature to probe SMBHBs and their coalescence: 
(1) precession of jet orientation and its acceleration in r adio galaxies during the in-spiraling 
of SMBHBs faegelman et all Il98pl : iLiu fc Chenll2007h . (2) optical periodic outbursts in 



blazars due to the interaction between SMBHB and accretion disk (ISillanpaa et al. 1 11988 



Liu et al. 


1995. 


2006; 


Liu k Wu 


2002: 



reorientation in X-shaped ra dio galaxies due to the exchange of angular momentum between 
SMBHB and accretion disk (ILiu 1 12004 ) . (4) two systems of broad emission lines (BELs) in 



at binary coalescen ce (ILiu et al. 



binary coalescence dMilosav 



quasars ( Boroson fc Lauer 2009i 



, (5) intermittent activity in double-double radio galaxies 



200 3J), (6) X-ray, UV, optical, and IR afterglow following 



2008: Schnittman &: Krolik 



emission lines ( Merritt et al 



jevic fc Phinney 1 120051 ; IShields &; Bonning 1 120081 : iLippai et al. . 
2008h , and (7) systematica lly shifted BELs relative to narrow 



(AGNs; iMadau fc Quataert 



2006 



2001 



Komossa et al.ll2008l) and off-center active galactic nuclei 



Loeb 1 120071 ) because of SMBH GW radiation recoil. 



All the above observational signatures require gas accretion disks around SMBHBs. 
The SMBHBs in gas-poor galactic nuclei are difficult to detect, because the SMBHs are 
dormant. However, a dormant SMBH could be temporaril y activated b y tidally disrupting a 
star passing by and accreting th e disrupted s tellar debris (jHills Ill975l ). A tidal flare decays 



typically as a power law t 5//3 (IRees 



several non-active galaxies (IKomossa fc Bade 



1988 



hmncy 



1999 



1989), which has been observed in 



Komossa 1 120021 : lHalpern et al. 1 12004 
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Esauei et al.ll2008l : iGezari et all 12008 . |2009|). Recently, IChen et all (120081 ) and IChen et al. 



( 120091 ) calculated the tidal disruption rate in SMBHB systems at different evolutionary 
stages, and found that it is significantly different from the typical rate for single SMBHs by 
orders of magnitude. This great difference of the flarin g rates enables o ne to statistically 
constrain the SMBHB population in normal galaxies (IChen et al.l 120081 ) . but identifying 
SMBHBs individually is still difficult at present. In this Letter, we investigate the influence 
of a SMBHB on the accretion of tidally disrupted stellar plasma and on the tidal flare. 
We show that the accretion of the stellar debris is interrupted by the SMBHB and this 
interruption can be taken as a key distinguishable observational signature for SMBHBs in 
gas-poor galactic nuclei. 



2. Tidal Disruption in a SMBHB system 



A star with mass m* and radius r* will be tidally disrupted by a SMBH with mass M, 
if it approaches the black hole (BH) within the tidal radius 



/ir,(M./m,) 1/3 



(jHills Ill975l ; iRees Ill988l ; iPhinney Ill989l ) , where /x is a dimensionless parameter of order unity. 
We focus on tidal disruptions of solar type stars with solar radius r Q and solar mass Mm by 



SMBHs with M, < 10 8 Mm, which are de t ected in the tidal-flare su rveys (IKomossa fc Bade 



1999: Komossa 



20021 ; lEsquei et al.M2008l : IGezari et al.ll2008 



2009). For more massive BHs, 



the tidal radii become smaller than the event horizon (jlvanov fc Chernyakovall2006l ). After 
tidal disruption at the tidal radius, the specific energy E across the star ranges from —E^ — 
AE to — Eb + AE, where Eb is the orbital binding energy of the star, AE = kGM,r^/r^ is 
the spread in specific energy across the stellar radius, and k = 1 if the tidal sp in-up of the 
star is negligible or k = 3 if the st ar is spun-up to the break-up angular velocity (IRees 111988 



Lacv et al. Ill982t iLi et al. 1 12002 ) . Since AE is (M./M*) 1 / 3 times larger than the internal 
binding energy of th e star and the interaction between the stellar debris is negligible after 
the tidal disruption (jEvans Sz Kochaneklll989l ). each fluid element moves ballistically. 



If we neglect the influence of the companion SMBH, the bound debris with E < should 
move in Keplerian orbits with eccentricity of about unity and fall back to the tidal radius 
after one Keplerian period, T = 2irGM,/(— 2E) 3 ^ 2 . The debris continuously returning to 
the tidal radius generates a falling-back rate m = (dm/dE)(dE/dt). For Keplerian orbits, 
dE/dt = dE/dT oc t" 5 / 3 . If we assume d m/dE following a constant distribution in the energy 
range — Eb — AE to —Eb + AE (see e.g. IRees Ill988l ; lEvans fc Kochanekj|l989l ; lLodato et al. 
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20091 ). the falling-back rate evolves as 



■ 

' > 1 mil 



t-T, 



D 



T ■ 

mm 



-5/3 



(2) 



for t > Td + T min , where To is the disruption time and 

2nGM. (2AE)~ 3/2 
0.35 5yr ^(^) 3/2 (^ 



- 1 mm — 



-1 



( M. 

VlO 7 M 



1/2 



(3) 



is the returning time for the most bound debris (typically <C AE). It is assumed that 
once the bound material comes back to the tidal disruption radius, it loses kinetic energy 
due to strong shocks because of the interaction between fluid elements and is circularized on 
a timescale muc h shorter than T to form an orbiting torus at a radius ab out r c = 2r t around 
the SMBH (e.g." 



Rees 


1988; 


Phinnev 


1989; 


Ulmer 


1999: 


Li et al. 


2002) 



dissipation of the shock energy and the accretion of the gas torus onto the BH will give rise 
to an X-ray/UV flare, decaying with (t — to) -5 / 3 - 

For a SMBHB system, the stellar-disrupting BH could be either the primary (with mass 
Mi) or the secondary (with mass M 2 ). Since the probab ility of stellar disruptio n by the 
secondary is relatively low if the two BHs are very unequal Jchen et al.ll2008l . l2009h . here we 
assume the primary BH to be the stellar-disrupting one. Because the interaction between the 
fluid elements is negligible before they come back to the tidal disruption radius, the SMBHB 
and each bound fluid element constitute a restricted three-body system. When the orbit of a 
bound fluid element is inside the SMBHB orbit, the system is called S-type. An S-type three- 
body system could be stable for a long time only when the system is hierarchical, namely, 
the system consists of an inner binary (the bound fluid element and the primary BH) on a 
nearly Keplerian orbit with semi- major axis Of and eccentricity of near unit, and an outer 
binary with eccentricity e and semi-major axis 3> Of in which the secondary BH orbits the 
mass center of the inner binary. The orbital change of a stable S-type system is negligible 
on the fluid-element dynamical timescale, so the falling-back stellar debris may interact with 
one another at the tidal disruption radius and finally be accreted with accretion rate given 
in Equation ([5]). However, for those less-bound fluid elements with af larger than a critical 
radius a cr , the orbit of the triple system becomes chaotic. For a triple system with an angle 
9 between the angular momenta of the inner and outer binaries, a cr can be semiempirically 
given by 



a b /a c 



2.8(l + g) 2/5 (l+e) 2/5 

x (l-er 6/5 (l-0.3#/180 c 



(4) 
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f lMardling Sz Aarseth Il200ll ). where q = M2/M1 < 1 is the BH mass ratio. 



Because of the nonlinear overlap of the multiple resonances flMardling 1 120071 ) in the 
chaotic triple systems, the fluid elements significantly exchange angular momentum and 
energy with the SMBHB and change their orbit dramatically on the dynamical timescale of 
the triple system, therefore would not return to the tidal radius to fuel the accreting torus and 
to form continuous accretion. Although a fraction of the fluid elements with chaotic orbits 
may return to the tidal radius on a timescale much longer than its Keplerian period, others 
would escape the system during the three-body interactions (see our numerical simulations in 
Section [3]). Therefore, Equation (J3J) implies that the continuous accretion of tidally disrupted 
stellar plasma stops with the accretion of the last fluid element with a < a cr , this occurring 
at a time 



-1/10 



-3/5 



x(l 



e) 9/5 (1 - 0.30/180°) 



-3/2 



(5) 



where Th is the SMBHB orbit al period. The eccentricity of hard SMBHBs is moderate 
( IMilosavljevic fc Merritt I l200ll ) and min or mergers with q -C 1 are the most common in 
the hierarchical galaxy formation model fjVolonteri et al.ll2003l ). We re- write Equation (jSJ) 
as T tr = 77Tb, where 77 is in the range (0.21,0.36) for e ~ and q <C 1, and independent of 
Ov For a SMBHB with semi-major axis = dh//3 residing in a stellar cluster with velocity 
dispersion a*, where 

GM,M 2 

(6) 



a h 



4(Mi + M 2 )al 

is the hard radius and (3 is the hardness, the accretion is interrupted at 



T tr = V T h ~ 6.2yra^ M 7 q 3 J?(l + q)~ 2 ( 



V 



0.25 



-3/2 



(7) 



where a 110 = aJllO km s~\ M 7 = Mx/10 7 M Q and g_ x = q/0.1. 



3. Numerical simulations and results 

Because the fluid elements of a tidally disrupted star move like test particles before 
falling back to the tidal radius, their evolution could be correctly simulated with scattering 
experiments. We use = 10 6 particles to resolve the stellar debris of a disrupted star and 
integrate their trajectories in restricted three-body systems. The particles are logarithmically 
sampled in the binding-energy range [E thd , AE], where E thd = G(M 1 +M 2 )/2.5ab is a binding 
energy corresponding to an orbital period of 4T), and AE = GMir^/r 2 is the spread in 
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specific energy for a non-spinning star. For a typical tidal- disrupt ion event with \Ej,\ <C AE, 
the ith particle has a binding energy Ei = AE(E thd /AEY i - Q -^' N , pericenter velocity V { = 
(2Ei + 2GM 1 /r t ) 1 / 2 , and mass m* = 0.5m*AEi/AE, where AE { = Eiln(0.5E thd /AE)/N is 
the size of the ith energy bin. 

The evolution of the three-body systems is computed in a frame centered on the mass 
center of the SMBHB with X—Y plane aligned with the SMBHB orbital plane. We assume 
e = in this work for simplicity, and that the secondary BH initially lies on the positive X 
axis and moves in the direction of positive Y . In each experiment, a test particle starts with 
velocity Vi from the pericenter at tidal radius r t about the primary BH, and the orbit is de- 
termined by three initial parameters (determined by the orbital parameters of the disrupted 
star at tidal radius) (1) the inclination angle 9 between the orbital planes of the particle and 
the SMBHB, (2) the longitude of ascending node Q, and (3) the argument of pericenter uj. 
With the initial parameters, we integrate the equation s of motion in the pseudo-Newtonian 



potential = GM,/(r — r g ) (jPaczyhski fc Wiita Ill980l ). w here r g is the Schw arzschild radius 



of BHs, using an explicit Runge-Kutta method of order 8 (lHairer et al.lll987l ). Stellar debris 
which falls back to within 2rt from the primary BH is assumed to be circularized via shocks 
and be accreted instantaneously. Therefore, if a particle with id i passes by the primary 
BH within a distance r acc = 2r t 0, the integration is stopped and the time tj is recorded. 
Otherwise, the integration continues until time 4T&. The accretion rate of the stellar debris 
is calculated using the recorded m; and U. 

Figured] shows the accretion rate of stellar debris in our fiducial simulations with q — 0.1, 
M 7 = 1, /3 = 10, (Xno = 0.9545, r* = r , m* = M Q , Q = 90°, and lu = 90°. Our results 
suggest that the interruption of accretion occurs at (0.15 — 0.5)T&, depending on 6. In 
Figure [lj a better agreement between the numerical value t tT and the analytical estimate 
T tr for 6 = 0° and 180° is by coincidence, because rj not only depends on 9 but also on Q 
and u. To illustrate this, we did 100 numerical simulations with random cos(0), f2, and to. 
The results show that r] ranges from 0.15 to 0.5 with a mean value 0.25, consistent with 
the analytical mean value 0.27 given by Equation (j3J). Figure [1] also shows that after being 
interrupted for about (1 — r))Tb, accretion recurs and "accretion islands" emerge discretely. 
The accretion rate at the islands is variable and can be larger than the corresponding value 
for single BH. The duration of interruption and the periods of the accretion islands decrease 
with 9, the shortest occurring at 9 = 180°. 

To investigate the effect of different on the results, we ran a scattering experiment for 
(3=1 and 9 = 0, keeping the other parameters as in Figure [TJ The result is given in Figure [2J 



1 We test our simulations with r aC c = l-5r{, 3rt, and 6rt and the results are nearly independent of r a 
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0.0 0.5 



1.0 
t/T b 



1.5 



2.0 



Fig. 1. — Accretion rates of gaseous debris in units of solar mass per year vs. time. Time 
starts at tidal disruption and is in unit of the binary orbital period Tb = 22.15 yr. The solid 
lines are the simulation results and the dashed lines give the standard ~ i~ 5 / 3 in single BH 
systems. T tr marks the interruption time estimated with Equation (jSJ). The upper, middle, 
and bottom panels are, respectively, for 9 = 0°, 90°, and 180°. 
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Fig. 2. — Accretion rates of gaseous debris vs. time for different at 6 — 0. The red and 
green solid lines correspond to /3 = 10 and /3 = 1, respectively. The other parameters are 
the same as those in Figure [U 
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Figure [2] shows that the interruption occurs at a slightly earlier time t tT /Tb as hardness (3 
decreases, which seems inconsistent with the analytical prediction given by Equation (JSJ). 
However, we should note that Equation (j5J) is obtained statistically by averaging over different 
Q and uj. When we did the simulations with different Q and uj, the averaged interruption 
time for different j3 resides in the predicted range. Accordingly, in the scattering experiments 
with different BH mass M\ (M7 = 0.1) and mass ratios q (q — 1, 0.02), we find that 77 weakly 
depends on q and M\. Our numerical results are consistent with the analytical estimate 
given by Equation (J5J), implying that the dependence of 77 on /?, q, and Mi is much weaker 
than on 9, Q, and uj. 



4. Discussions 

We investigated the accretion of tidally disrupted stellar debris in SMBHB systems. For 
simplicity, we assumed for the tidal debris gas: (1) a constant distribution of mass in binding 
energy at the tidal radius, (2) ballistic motion and negligible interaction of the fluid elements 
after tidal disruption, and (3) instant circularization of the debris gas probably due to shocks 
because of interaction between the fluid elements returned to within a radius two times the 
tidal radius. T he first two assumptions have been justified by numerical hydrodynamic 



simulations (e.g. lEvans fc Kochanekl Il989l ; lLodato et al.ll2009i ). but the third one needs to 



be verified by hydrodynamic simulations capable of capturing strong shocks. With these 
three assumptions, the fluid elements and SMBHB compose restricted three-body systems, 
so we can investigate the accretion of the gaseous stellar debris by analy zing the stabilit y 



of the three-body systems using the resonance overlap stability criterion (IMardling 1 120071 ). 
Because of the chaotic nature induced by the nonlinear overlap of several orbital mean motion 
resonances, the fluid elements inside the chaotic regions significantly change their orbits on 
a dynamical timescale and do not return to the tidal radius to fuel the BH, leading to the 
interruption of accretion. We also investigated the evolution of the stellar debris using three- 
body scattering experiments. Our results obtained both analytically and numerically show 
that the accretion rate of the debris gas decreases with a power law M ~ t _5//3 until a critical 
time T tr . At time t > T tr , the accretion pauses until about one SMBHB orbital period T h . 
T tr relates to T\, with T tr = 77Tb, where 77 is typically 0.25 and in the range 0.15 < 77 < 0.5, 
depending on the initial orbital parameters (9, Q, uj) of the fluid elements at tidal radius and 
on the eccentricity of the SMBHB, but being nearly independent of the SMBHB semimajor 
ab, BH masses, and mass ratio q. This suspension of accretion would result in an interruption 
of the tidal flare, although the residual accretion disk may still radiate weak optical, UV, 
and X-ray emission during the interruption. 
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Our numerical results indicate that the accretion of the gaseous stellar debris restarts 
at a time T r , leading to a flicker of flare. The exact recurring time, T r = £Tb, depends on the 
initial orbital parameters 9, Q, and uj, but our numerical results suggest that £ is of order 
unity and 1 < £ < 1.5. The interruption timescale of the tidal flare in a SMBHB system is 

AT = T r — T tr ~ (£ - 77)T b . (8) 

Our numerical simulations suggest that the accreted plasma during the discrete accretion 
consists of both the tidal debris gas falling back to the tidal radius after one Keplerian time 
and a fraction of those fluid elements with chaotic orbits falling back to the accreting torus on 
a timescale longer than its Keplerian time. When T tr and AT is determined observationally, 
constraints on 77 and Tb could be made if we take £ — 1. 

In our simulations, we assume that the orbital binding energy Ef, of the tidally disrupted 
star is negligible compared to the spread in specific energy AE. For such tidal-disruption 
events, Equations ([3]) and ([7]) imply that the standard power-law decay and the interruption 
of tidal flares are detectable if 



/?<Anax « lSOa^M^V^l + g)- 4 / 3 



k 

2 

2/3 



T) \ 2 /3 (R*\ f m* , 

d257 It:J ImIJ • (9) 



If the tidally disrupted star is initially y ery bound to one of the binary BHs so that \E t 



6 



is comparable to AE f lChen et al.1 120091 ) . the time T mm would be much shorter and the 



interruption of tidal flare is detectable even in an ultra-hard SMBHB with (5 3> /3 max . For a 
SMBHB with /3 ~ 100, the semi-major axis of the binary ab is about ~ 9.3 x 10 2 r g g_i(l + 
q)~ 1 aii Q . The GW radiation emitted by such SMBHBs could be detected by PTA. 

An interrupted tidal flare could be caught if T tr is shorter than the mission duration of an 



transient survey t sv ~ 1 — lOyr, which corresponds to SMBHBs with j3 ~ (5 — 100) 



Here we use the M.-a* relation ( Tremaine et al.ll2002h . Because a hard SM BHB with q = 



1 



spends most of its life time at (3 ~ 5 — 100 f Yu II2OO2I ; ISesana et al. 112008 ). we have a good 



chance to detect them with upcoming transient surveys. However, if one wants to catch an 
interrupted tidal flare from SMBHBs emitting strong GW radiation (f3 ~ 1000), the time 
resolution of the survey should be < 0.1 yr. Because a SMBHB spends a small fraction of its 
lifetime at this stage, high-sensitivity and deep transient surveys are needed to accumulate 
many more tidal-disruption events. 
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